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Abstract 

We consider multiple sums and multi-integrals as tau functions of the BKP hierarchy using 
neutral fermions as the simplest tool for deriving these. The sums are over projective Schur 
functions Qa for strict partitions a. We consider two types of such sums: weighted sums of Qa 
^ ^ over strict partitions a and sums over products QaQ-y ■ In this way we obtain discrete analogues 

^ ■ of the beta-ensembles (/3 = 1,2,4). Continuous versions are represented as multiple integrals. 



Such sums and integrals are of interest in a number of problems in mathematics and physics. 



1 Introduction 

This work deals with certain multiple sums ( fl2.ip -( l2Tl) . f l2.20p - fl2.23p ) and multiple in- 
tegrals (f l3.ip -f l374|) . fl3.20p ) that are of interest in a number of problems in mathematics 
and physics. The multiple sums appear in models of random partitions (as developed by 
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A.Vershik's school; see [T] for a review), and random motion of particles [2] (see also |1] 
for a review). The multiple integrals yields certain analogues of /3 = 1,2,4 ensembles of 
random matrices, and also an analogue of the two-matrix models We consider defor- 
mations of the measure defined in terms of four semi-infinite sets of parameters and relate 
the generating function Z to the coupled two-component BKP hierarchy. The main tool 
is the use of 1- and 2-component neutral fermions, which provide links with integrable sys- 
tems known as neutral BKP hierarchies [3] , [0] . One version of these was first introduced 
in [3], [0]. In this work, we use another version introduced in [7]. 



2 Sums over projective Schur functions 

In the following, we consider sums over strict partitions , which will be denoted by Greek 
letters a, /3. Recall |8] that a strict partition a is a set of integers (parts) (ai, . . . , a^) 
with «! > ■ ■ ■ > ftfc > 0. The length of a partition a, denoted i[pi)^ is the number of 
nonvanishing parts, thus it is either k or k — \. Let DP be the set of strict partitions 
(i.e., with distinct parts). We also need a subset of DP, which will be denoted DP^ and 
consists of all partitions of the form (ai, a\ — 1, 03, 03 — 1, ... , Ofc-i, — !)• Consider 
the following sums (for L G N+, t := (^1,^3, . . . ), t* := {t\,t\, . . . ), t := i^xM, •••))• 

5o(t,L) := (2.1 

oigDP 

5i(t,t*) := e-'^-^^'^Qadt) (2.2 

aGDP 

^2(t,t,r) := e-''"^**^Qa(|t)Q„(it) (2.3 

qGDP 

5oo(t,t,L) := Y (2-4 
53(t,A^) := Y (2-5 

aSDP 

S,{tX) ■■= Y e-^'-^^^^Qadt) (2.6 
S,{t,t,D) := Y QX2i)Da,pQpC2^) (2.7 



l(c.)=l{fi) 



Here, the projective Schur functions Qa are weighted polynomials in the variables 
t = (ti, t3, ^5, . . . ), degtm = m, labeled by strict partitions (See [8j for their detailed 



2 



definition.) Each Qa{\'t) is known to be a BKP tau function [5], [6]. (Tliis was a nice 
observation of [H], [ID])- Tlie fact that only odd subscripts appear in the BKP higher 
times t2m-i is related to the reduction from the KP hierarchy. The coefficients Ua are 
defined as 

k 

U^:=Y,Ua,. (2.8) 

i=l 

where 

[/„:=[/W_ n'"4-lnn!, n G N+ (2.9) 

rriT^OjOdd 

for some given set of constants {Un^}. 

The coefficients on the right hand side of fl2.5p are determined in terms a pair 
{A, a) where A is an infinite skew symmetric matrix and a an infinite vector. For a strict 
partition a = (ai, . . . , a^), where ak > 0, the numbers A"^ are defined as the Pfaffian of 
an antisymmetric 2n x 2n matrix A as follows: 

Al := Pf[i] (2.10) 

where for k = 2n even 

= -iji := l<i<i<2n (2.11) 

and for A; = 2n — 1 odd 

7 7 iAaa if 1 <«<?■< 2ra-l 

Ai = -Aji := < " ' ~ ~ (2.12) 

|^a„, if 1 < z < j = 2n. 

In addition we set Aq = 1. 

The coefficients Da,i3 in (12. 7p are defined as determinants: 

D,,^ = det (D^^^^^) (2.13) 

where D is a given constant infinite matrix. 

Remark 2.1. Series (j2.ip . (j2.2p . (j2.3p and (j2.4p may be obtained via specializations of A'^ in 
the series ()2.5p . If we put 

Am = -e"^'"~^"sgn(n - m), = e"^" (2.14) 

we obtain (j2.2p . If we further choose [/„ = and C/„ = +oo for n < L, n > L respectively, we 
obtain (|2.ip . If we set 

Anm = 2^ ^"Q(n,m)(f t)5 O-n = ^ '^"Q{n) (^"t) (2-15) 



we obtain (j2.3p . Choosing again Un = 0, [/„ = +00 for n < L,n > L, we obtain (j2.4p . The 
series (j2.6p is obtained from ()2.5p by taking 

The sums (I2.4p and (I2.3p may be also obtained as particular cases of (12. 7|) by putting 
to get ([231). 

All these sums are particular examples of BKP tau functions, as introduced in [5], 
defining solutions to what was called the neutral BKP hierarchy in 0. They may be 
further specialized if we choose t = too:= (1,0,0,...). Then 

k 

Q^{lt^) = A*{a)f[^, « = («!,..., afc) (2.18) 

1=1 

where 

A*(a)=A*(a):= J] (2-19) 

0<i<j<k ■' 

For this specialization, we have 

Si{t^X) = J]A*(«)n — (2.20) 



ail 

leDP i=l 



52(too,too,r) = ^A*(a)2n^l^ (2.21) 



^ g-C/,,{t*)-C/<., + i(t*) 



oGDP' *=1 ^ ^ 

55(too,too,t*) = 5^ - 5^ n ^ (2-23) 

«(a)=f(/3) = fc 

Here DP' is the set of all strict partitions (ai, q;2, • • • , «Af > 0) with the property 
ai > CKj+i + 1, z = 1, . . . , — 1, and 

i<j<N + "j') + "^'^ ~ ^) 

If we replace the term A* in the above by the Vandermonde determinant A, the 
resulting sums fl2.20p - fl2.23p may be viewed as discrete analogues of matrix models (see 

my 
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Applications. The following are some examples of applications of the above sums. 

(I) The sum (12.41) was considered first by Tracy and Widom in [12] in a study of the 
shifted Schur measure. 

(II) Sums (12.21) and (12. 3p may be viewed as generalizations of hypergeometric functions for 
the case of many variables. For example, a generalization of the hypergeometric function 
of type pF^ may be obtained from expressions (12.21) and (12. 3p by defining the parameters 
Un in terms of Gamma functions as follows: 

U., = log "r'^i":?'! (2.25) 
11^=1 r(n + bi) 

for some set of p + r constants {ttj, bi}. Sums of type (12.31) were considered in [13], while 
sums (12. 2 p are new. (This case will be considered in detail elsewhere [2Z]-) It may be 
shown that both series (12.20 and (12. 3p may be expressed as Pfaffians of matrices whose 
entries are expressed via pF^ and share many properties with the usual hypergeometric 
functions pF^. Analogues of basic hypergeometric functions may be obtained in a similar 
way. 

(III) Consider models of random strict partitions a, where the relative weight Wa is given 
by one of the following: 

(A) 

W„ = A^g„(|t), (2.26) 
where A'^ = {A, a) and t = {ti,ts, . . .) are parameters of the model 
(B) 

= e-^-^^^^^Q^dt), (2.27) 
where the parameters are t = (ti, ^3, . . . ) and U = {Uq, Ui, . . .). 
(C) 

= e-^'^(**)Q«(it)Q,(it), (2.28) 

where t = (ti, ^3, . . . ), t = (ti, ^3, . . . ) and U = {Uq, Ui, . . .) are independent parameters. 
(Note that models (B) and (C) are particular cases of model (A).) 

The series S3, Si and may then be be viewed as normalization factors (partition 
functions) respectively for models (A),(B), and (C). Similarly, series 5*5 is a partition 
function for a model of strict bi-partitions. 

(IV) Series Si and S2 were used in [18], where random oscillating Young diagrams related 
to strict partitions were considered. 

Remark 2.2. The fermionic representation for these models allows their correlation functions 
to be computed in standard ways. (See e.g. refs. [H], [H], [16] and [TT].) 
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Neutral fermions. To construct tau functions of the BKP and two-component BKP 
hierarchies (see [6], [?]) we need neutral free fermions {(f)f \ (f)f^}i(zz, satisfying the anti- 
commutation relations 

= i-irSn,-mSa,b, a, 6 = 1, 2. (2.29) 

For one-component BKP only the first component 0„ := (pn \ r?, G Z is used. The fermionic 
Fock space will be chosen as in [7j) (as opposed to [6]). Namely, the action of neutral 
fermions on vacuum states is defined by 

(^i")|0) = 0, (O|0Lt = O, n<0, (2.30) 

4")|0) = -^|0), (014'^) = -^(0| (2.31) 

where relation (12.311) is chosen as in [7|). (See also the Appendix A in the arXiv version 
of ^18i| for some details on links between [6] and [7].) 
For linear combinations 



Wk 



EE^n^i'^^' k-1,2,..., (2.32) 



Wick's Theorem implies, for arbitrary such products of an even number of WkS 

{Wi ■ ■■W2n) = E Sign(0-)(«;^(l)t(;^(2)) ■ ■ ■ («^a(2n-l)«^a(2n)) =: Pf (( (^iWj ) ) j<2n 

where the sum is over the permutation group S2n- 
We also need the Fermi fields 



(2.33) 



0(«)(z) = E0i"^^"> « = 1,2 (2.34) 



To evaluate integrals we use 



(O|0(^)(zi)0('^)(z2)|O) = i (^^-^ ) 6a,. (2.35) 



For m G Wick's theorem implies 



(O|0('^)(;.i)0^'^n^2)---0('^n^J|O) ={l)' Aliz) (2.36) 



2 

Note that, because of (12.311) . this expectation value is nonvanishing for m odd. 
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BKP tau functions. The general tau function of tlie two-component 2-BKP liierarcliy 
is expressed in fermionic form as 

= (0| r«(t(i))r(2)(t(2)) eSa..E„.„ A^^U'nUl^' f «(t(i))f (2)(t(2)) |0). (2.37) 

where 

r('^)(t('^)):=exp Yl Kt^n^ r('^)(t('^)):=exp 5^14") (2.38) 

n>l, odd J^^l, odd 

■= \Y.(-^y^'<^^^^-ln- (2-39) 

The array of numbers A = {A^^, a, 6 = 1,2; n,m G Z} is the data that determine 
the two-component 2-BKP tau function. The four sets of independent parameters, t^"^ = 
{t^i\t^^\ ...), t^"^ = {t^i\ ■■■), ior a = 1,2, are called higher times of the hierarchy. If we 
fix any three sets, the fourth will be the higher times of the usual BKP hierarchy [7]. In this 
sense (12.371) may be viewed as a four coupled BKP tau function. If we set t^^^^ = t*^^) = 
we obtain the two-component BKP hierarchy, as described in [7j. 

For the one-component case we omit the second component and the superscripts. The 
2-BKP tau function is then 

T^^^^it, t) = (0| r(t) e^-.'-ez A,,^<j>„<p^ Y{t) |0) (2.40) 

and the usual neutral BKP tau function may be written 

r^^-P(t) = (0| r(t) e^".™« |0) (2.41) 

A remarkable example of a BKP tau function was found in [9] ; namely, Schur's Q-functions 
Qa themselves |8j. These may be expressed fermionically as 

(0| r(t) 0^,0,, ■ ■ ■ |0) = 2-i^(")Q,(|t) (2.42) 

where ai > a2 > ■ ■ ■ > > 0. The set (ai, . . . , a2Af) and the partition a are related 
as follows: in case > 0, a = (ai, . . . , q;27v) and £(a) = 2N , while in case = 0, 
a = («!, . . . , a2N-i) and i{a) = 2N — 1. 

Fermionic representation for sums. The formulae below show that sums fl2.ip - fl2.6p 
are BKP tau functions. First, we have 

S3{t,A^) = J2 ^a<5a(|t) = (O|r(t)e2S">™>o^™'^"'^'" + 2E„>o«n0n</'o |o). (2.43) 

aeDP 
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In view of Remark 12. II we have fermionic representations for (12.11) . (12. 2p . (12. 3p . (12. 4p . In 
particular 

So{t,N)= ^ g„(it) = (0|r(t)e2^^>">'">«'^"'^™ |0) (2.44) 

aeDP 
ai<N 

Si{t,t*) = e"^"^**^<5«(|t) = (0|r(t)T(t*)e2S->->o<^"'^™ |0) (2.45) 

S^{t, t*) = J2 e"^"^**^(5a(|t) = (0| r(t) T(t*) e2^">o ^"'''"+' |0) (2.46) 



where 



T(t*) := exp l-J2Unit*)(pn(p^n] (2.47) 

\ n>0 / 

Evaluating at t = too gives a "solitonic" representation for Si{t = too,t*) : 
5i(too,t*) = i(0| r(t;) e2E„>™>o e-^l°'-^i°V(n)</'M+2E„>o<^-^"'"</'(n)<^o r(t*_) |0) (2.48) 

where Fermi fields (12.341) are used (see [H]). This shows that this series is a 2-BKP 
tau function with respect to the higher times defined as t'j_ := {tl,tl, . . . ) and t*_ := 
{—t*_^, —t*_^, . . .). Here c is a normalization factor given by 

(o|r(t;)r(r )|o) = b{t*_,ti). (2.49) 

Remark 2.3. Another fermionic representation for (I2.3p .( f2^ was obtained in |13] . 
Finally, we have 

S,{t^'\t^'\,D) = (0|r«(t«)r(2)(t(2))e2E„,„>oA.^</'l^V^^^|0). (2.50) 

3 Multiple integrals 

Let diy be a measure supported on a contour 7 on the complex plane. We take 7 as either 
of the following two contours: 

(A) An interval on the real axes < z < 00 

(B) A segment of the unit circle: given hj z = e^'^ , < < 9, < 9 < n, 
Consider the following A^-fold integrals: 

h{N):= I ■■■ I \/\*{z)\\{dv{z,) (3.1) 
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N 



■■ / \A*{z)\'l[du{z,) (3.2) 
h{N) := /"••■/" A*(z)a'=(z) II dv{zi) (3.3) 
h{N):= I ... I \/\\z)\'\{dv{z,) (3.4) 



where, as before, 



N 



A*(.) = n 



Zi Zj 



Zi -\- Zj 

1>J -I 



The notation a^(z) is analogous to f l2.10p . denoting the Pfaffian of an antisymmetric 
matrix a: 

a'=(z) := Pf[a] (3.5) 

whose entries are defined, depending on the parity of A^, in terms of a skew symmetric 
kernel a{z, w) (possibly, a distribution) and a function (or distribution) a{z) as follows: 
For N = 2n even 



a. 



-dji := a{zi, Zj), 1 < i < j < 2n (3.6) 



For = 2n - 1 odd 



{a{zi, Zj) if l<i<j<2n— 1 
(3-') 
a{zi) if 1 < i < j = 2n 

In addition we define Oq = 1. 

Integrals /i,/2 and J4 may be considered as analogues of /3 = 1,2,4 ensembles [3]. 
They may be obtained as particular cases of I3 as follows: 
Integral h{N) is a particular case of IsiN) where in the (A) case 

a{zi, Zj) = sgn{zi - Zj), a{z) = 1 (3.8) 

while in case (B) 

a{zk, Zj) = e~ 2 sgn(v9fc - (pj), a{z) = e" 4 , (3.9) 
with ifi = argZj. To prove this we use: 
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Lemma 3.1. 

Pf [sgn{zk - Zj)] = sgn A*(z), Zk E 

Pf [sgn(y., - ^j)] = sgn (^e-f (^^-^)A*(z)^ , 

where k, j = 1, . . . , N . 

Integral hiN) is obtained from l3{N) by setting 



Zk = e 



(3.10) 
(3.11) 



Zi Zj / \ 1 

a{z) = 1. 

Zi ~\~ Zj 



We use the fact that 



A*{z) = Pf 



Zj 



^Z{ -\- Zj _ 

Integral I^lN) is obtained from l3{2N) as follows. In case (A) we set 



a{zi, Zj) 



and in case (B) we set 



^(^25 Zj) 



d_ 



■5{ip, - tpj). 



(3.12) 
(3.13) 

(3.14) 
(3.15) 



The integrals containing A* may be compared with those defining the partition func- 
tion of the so-called supersymmetric matrix integrals, see [19]. Integral I2 defines the 
partition function of the so-called Aq model (see [20]), the Coulomb gas model with re- 
flection [21], the ID Ising model [22], and correlation functions in the 2D Ising model [23] . 

To relate these integrals to the 2-BKP hierarchy we introduce deformations Ii{N) — )■ 
/j(A^;t,t) through the following deformation of the measure 

(3.16) 



where 



and 



du{z) diy{z\t,t) = b{t, {z})b{-t, {z-^})du{z) 

V ^ n 

6(s,t) = exp -Sntn 



odd 



, , , 2z^ 2z^ 
{z}^i2z,-,-, 



(3.17) 



(3.18) 



Below, we show that the generating series obtained by Poissonization (the grand par- 
tition function) 



00 N 

Z,{fi ; t, t) = b{t, t) t, t) ^ , ^ = 1, 2, 3, 4, 



(3.19) 



N=0 
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are particular 2-BKP tau functions fl2.40p . 

We also consider the following 2A^-fold integrals: 

N 

/5(iV;tW,t(2),t«,t(^)) := / A^(z)A;.(y)nM^.y.|t«,t(2\t«,t(2^ (3.20) 

i=i 

where 

dK^,y|t«,t(2),t(i),t(2))= (3.21) 

b{t^'\ {z})b{-i^'\ {z-'mt^'\ {ym-i^'\ {y-'})du{z, y) 

(here di>{z,y) is an arbitrary bi- measure), and show that the generating series 
Z5(^;t(i),t(2),tW,t(2)) = 6(t(i),tW)6(t(2),t(2))^ J5(iV;tW,t(2),t(i),t(2)):^ (3.22) 

N=0 

is a particular case of the two-component 2-BKP tau function (12.371) . 
Remark 3.1. Note that 

Z2 (m ; t, t) = Z5(m ; t(i) , t(2) , t(i) , t(2)) (3.23) 

if 

du{z,y) = 5{z -y)diy{z)diy{y), t = tW+t(2), t = tW+t(2). (3.24) 

The integrals Zi(/i;t,t), Z2(/i;t,t), Z4(^;t,t) and Z5(/i ; t^, t^^), t^, t^^)) may be 
obtained as continuous limits of 5*1 (too, t*), 'S'2(too5 too, t*), 'S'4(too,t*) and 5*5(100, too, t*), 
respectively. 

Fermionic representation of the integrals. To obtain the fermionic representation 
for the integrals above we apply ( 12.36^ . Expanding the exponentials and applying Wick's 
theorem to each term in the sum gives 

Z3(/i;t,t) = (0| r(t) ^1^^ ^'^i^,y)'t'i^)Hy)duiz)du(y) ^2t,J^aiz)Hz)^odu{z) |q^^ ^-3^25) 

Note that 

ZsiO ; t, t) = 6(t, t) = (0| r(t) f (t) |0) (3.26) 
was written above in (13.171) . In particular we obtain 

Zi(/i-t t) = (o| r(t) e'^^'^^ I-y fj'^s^^'^^''')~'^^yy')'f'^''')'f'^y^'^'^^^^'^''^y^ e'^>^^ |o) (3 27) 
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where q = 1 and c; = '^{z) , 2; G 7 is a parameter on 7, which is equal to q{z) = z when the 
integration contour is IR+, while in the case 2; = e**^ G 7, < (/? < ^, we set (;{z) := (f and 
q = e 4 . 

Similarly, we have 

Z2(^;t,t) = (0|r(t)e^'^^ASf<^W<^(3/)'^'^(-)'^''(2^)e2^A'^(^)'^o^^^^^ |0) (3.28) 
As a specialization of (13. 25 p . by choosing a{z,w) as in (13.141) and a{z) = 0, we obtain 

Z4(/it,t) = (0|r(t)e^^^-r"^'^(")'^"(")f(t) |0). (3.29) 
Finally, in terms of two-component fermions we have 
Z5(/i;t«,t(2),tW,t(2) 

= -(01 r«(tW)r(2)(t(2)) e^^^// ^^^HzH^'Hy)Mz,y) f {i)(t(i))f (2)(t(2)) |o), (3.30) 
c 

where c := ^5(0 ; t^^\t^'^\t^^\t^'^^) is the normalization factor. 

Formulae (^M> ,(^^ ,<^^ , and (K30\\ should be compared with (I233]),([225D,([226D, 
and (I2.50p . respectively. The fermionic representations of integrals Zi and Z4 are similar 
to the results of [24j for ensembles of real symmetric and self-dual quaternionic random 
matrices. 

4 Discussion 

We have presented five types of multiple sums and multiple integrals which are related to 
tau functions of BKP type hierarchies (BKP and two-component BKP tau functions for 
multiple sums, 2-BKP and two-component 2-BKP tau functions for multiple integrals). 
Certain of these sums and integrals have known applications in mathematics and physics; 
we believe that all of them may prove to be of use in various probabilistic models. The 
techniques of free fermion calculus and integrable systems may be applied to study various 
properties of these sums and integrals. Multicomponent BKP tau functions [7] may be 
further used to construct models of Pfaffian processes (cf. [26]). The results of this work 
should be be compared with analogous results for the so-called charged BKP case (see j^) 
studied in pi] and [27]. 

5 Appendix. An application of the integral Z2{N) 
(Harry Braden). 

Some of the considerations of this paper were motivated by the following observation of 
Harry Braden ( [23]). 
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"The apphcation of the integrals I have in mind deals with an Ising model corre- 
lation function when there is a thermal perturbation from critical temperature. In the 
scaling limit this is described by a system of Majorana fermions and makes connection 
with the fermionic representation under consideration. McCoy, Tracy and Wu showed this 
limit was described a massive (m = T — Tc) field theory whose correlations were governed 
by a radial Sinh-Gordon equation that under a change of parameters is Pm. Another 
approach to the same correlator is via form factors. This approach yields the correlator in 
terms of an infinite sum of multiple integrals. For the case at hand this gives a Euclidean 
correlation function like 

G(r) :=< O(x)O(0) > = f; / < 0|O(x)|/3i . . . >< . . ./3i|O(0)|0 > 

n!(27r)" 

= E / %4§im---/3n)rexp(-mrf;coshA). 
Here /3j are rapidities and x = {xq,xi), r = ^/xl + x\. If we use the minimal form factor 



Fr"(/3i.../3„) = ntanh 



and appropriate normalizations we get 



A - A- 



As I havent been specific about the precise correlator, we also are interested in 

n=0 i=l ^ ' i<i ^ ■' ' 



The identity 



\ 'T' ■ —I— T' ■ I /^'T'^ T" J-J-V'T'- — U ■ 

\ »AV 2 I »AV ^ / £^ tXj \ • • • ikj nn \ 1 I lAv •! 



provides a connection with a Fredholm determinant. I have looked at expansions of (15. ip . 
The difficulty is getting a convergent expansion. The first term is in terms of i^o(^)- Such 
an expansion will provide one for (a particular) solution of the Painleve equation." 
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